The Resonating- Valence-Bond Ground State of Li Nanoclusters 
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We have performed Diffusion Quantum Monte Carlo simulations of Li clusters showing that 
Resonating- Valence-Bond (RVB) pairing correlations between electrons provide a substantial con- 
tribution to the cohesive energy. The RVB effects are identified in terms of electron transfers from s- 
to p-like character, constituting a possible explanation for the breakdown of the Fermi liquid picture 
observed in recent high resolution Compton scattering experiments for bulk Li. 
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Lithium, and other alkali metals, have been modeled 
as free-electron like systems - an electron gas permeated 
by ions [ij. Nevertheless, in recent years, experimen- 
tal and theoretical investigations of Li have revealed a 
more complex phase diagram |^. Even at ambient 
pressure, the electron momentum density cannot be de- 
scribed adequately in terms of Fermi liquid theory, be- 
cause pronounced deviations have been observed in bulk 
Li in recent high resolution Compton scattering experi- 
ments 0, 0] ■ Such deviations from the standard metallic 
picture can be ascribed to the possible existence of signifi- 
cant pairing correlations in the ground state @| . Bonding 
properties have also revealed that Li behaves like a " bad" 
free-electron metal Q ■ 

The notion of stabilizing the metallic state through the 
creation of a resonant valence bond (RVB) state involv- 
ing the metallic orbitals dates back to the early works of 
Pauline^ who first applied this picture to the Li ground 
state [Sj. In 1987, Anderson proposed the RVB wave 
function as the natural ground state for the high tem- 
perature superconducting materials [9j , arguing that this 
picture is capable of capturin g rn any aspects of the phase 
diagram of the cuprates More recently, it has 

been shown fj] that Pauling's RVB idea cannot account 
for all of the properties of metals that depend on the 
existence of a Fermi Surface (FS). However, since the na- 
ture of the FS in bulk Li has been questioned by the high 
resolution Compton scattering studies , it is natural to 
ask if the RVB paradigm might provide a viable model 
of the Li ground state. 

This work demonstrates that RVB pairing correlations 
between electrons provide a contribution to the total en- 
ergy of Li clusters of about 20 meV or greater per atom. 
The pairing correlation effects modify the electron mo- 
mentum density distribution [i] and therefore provide a 
possible mechanism for the breakdown of the Fermi liquid 
picture in bulk Li. 

A number of authors have performed calculations on 
bulk Li as well as Li nanoclusters Q , but an implementa- 
tion of the RVB model for Li clusters utilizing QMC sim- 



ulations has not been attempted. A previous RVB study 
of small Li clusters [I3[ did not possess the accuracy of 
modern QMC methodologies. Here, we report QMC cal- 
culations of correlation effects beyond the limitations of 
the standard Jastrow-Slater wave function (WF) for the 
Fermi liquid ground state [l^ by employing the RVB. 
Specifically, we have used the Jastrow Antisymmetrized 
Geminal Product (Jastrow-|-AGP, or J AGP), developed 
by Sorella and coworkers [11, 0, |l3, H El , as a QMC 
variational ansatz for the RVB. We have performed Dif- 
fusion Monte Carlo (DMC) calculations to obtain precise 
estimates of the energy, this technique being limited in 
accuracy only by the nodal structure of the variational 
ansatz [20|, l2l| . For a useful DMC study of bulk Li, see 
[2^ 1. For all our Li clusters, we obtain a distinct nodal 
structure improvement of the cohesive energy in compar- 
ison to standard JS WFs. An eigenvalue analysis of the 
pairing wave function further confirms the RVB nature 
of the ground state. 

Our calculations employ the J AGP WF, defined 
as the product of a Jastrow term J and an anti- 
symmetrized determinant 'i'AGP- jagp{i~i, ■■■,fN) = 
^AGp{ri,--,rN)J{ri, ...,rN)- The determinant *agp 
is constructed as Iiagp = i[«'(r|, rj;) • • • $(rJy/2, ^-jv/a)], 
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where A is the antisymmetrization operator, <i>(r|,rj) is 
the pairing function, and N is the number of electrons 
in the system {N must be even for this formulation; see 
[2^ for an extension to open-shell systems). The pairing 
function $ is defined as <I>(r^,r^) = ■0(r^, r^)I/\/2(| ti 
) ~ I iT))i where the spin part is a singlet. The spatial 
part ijj - the geminal - is represented by a pairing expan- 
sion over a local single-particle basis set {fpi}, i.e. 
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Here, indices i and j run over different orbitals (covering 
all nuclear sites), which are expanded in a gaussian basis 
set centered on the nuclear positions The Jastrow 

factor J = J1J2J3 is composed of an electron- nuclear 
(Ji), an electron-electron (J2), and an electron-electron- 
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FIG. 1; (Color online) Cohesive energies per atom for Li„ 
with n = 2, 4 and 8, representing systems of dimension d = 1, 
2 and 3, respectively, calculated using the diffusion Quantum 
Monte Carlo all-electron method. The plot compares the RVB 
nodal structure against the standard IfF nodal structure. The 
error bars are comparable to the size of the points and lines 
connect points to guide the eye. The inset shows the difference 
in cohesive energies between the models. 



nuclear ( J3) term; it guarantees that the cusp conditions 
are satisfied and it allows the correct charge distribution 
in the system. The J3 term is constructed in a form 
similar to the pairing function of Eq. [1] ^16]. The Jas- 
trow parameters, A parameters, gaussian (Slater) orbital 
exponents, and orbital coefficients have been optimized 
by energy minimization using the methods described in 
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We have compared the JAGP wave function with three 
types of single-determinant JS wave functions: one in- 
volving a standard Slater determinant of Hartree-Fock 
(HF) orbitals (abbreviated JS-HF); one utilizing a pseu- 
dopotential to replace the core electrons and utilizing 
Kohn-Sham orbitals for the valence electrons calculated 
within the Density Functional Theory (DFT) using the 
Local Density Approximation (JS-LDA); and a wave 
function defined as the limiting case of the JAGP wave 
function in which the occupation of the virtual orbitals is 
forced to be zero (JS-VO). The latter wave function pos- 
sesses the standard JS form. The geometries of all our Li 
clusters were optimized using the software Jaguar with 
a 63 111 G basis set and the B3LYP exchange-correlation 
potential [26 



27| 



For the all-electron JS-HF wave function, we examined 
clusters of 2, 4, and 8 Li atoms modeled with HF or- 
bitals obtained from the software Jaguar [26| . The wave 
function included the Ji and J2 terms, which were opti- 
mized using an improved version of the Stochastic Gradi- 
ent Optimization (SGA) 28] method. Our correspond- 
ing DMC results are given in Table U for the dimer. 



N 


rpJS-HF 

^DMC 

(Hartree) 


rp.JAGP 
JUDAIC 

(Hartree) 


coh 

(eV) 


■pJAGP 

(eV) 


(meV) 


2 


-7.49593(8) 


-7.4971(1) 


0.491(2) 


0.522(3) 


32(5) 


4 


-7.5036(1) 


-7.50485(13) 


0.699(3) 


0.733(4) 


34(6) 


8 


-7.5130(1) 


-7.51440(23) 


0.955(3) 


0.993(6) 


38(9) 



TABLE I: DMC all-electron calculations comparing the JS- 
HF and the JAGP wave function for Li clusters. The first 
two columns present the total energy/atom. The next two 
columns show the cohesive energy/atom, and the final column 
is the difference between these two cohesive energies. 



N 


rpJS-LDA 
^DMC 

(Hartree) 


jT,JAGP 

(Hartree) 


r,JS~LDA 
^coh 

(eV) 


jT,JAGP 

(eV) 


AEcoh 

(meV) 


2 


-0.21318(4) 


-0.21534(8) 


0.45932(7) 


0.51625(9) 


56.9(2) 


4 


-0.22100(3) 


-0.22357(5) 


0.67212(6) 


0.74020(6) 


68.1(1) 


8 


-0.23152(3) 


-0.23225(1) 


0.95839(6) 


0.97640(2) 


18.01(8) 


20 


-0.237750(5) 


-0.23893(1) 


1.12791(4) 


1.15817(2) 


30.26(6) 



TABLE H: Same as Table HI except here DMC pseudopoten- 
tial calculations are used to compare the JS-LDA and the 
JAGP wave function. This data set includes Li2o . 



the planar cluster Li4 and the three- dimensional Lig. 
For all three clusters, we obtain distinct corrections of 
about 30 meV/atom for the cohesive energy. Fig. [T] il- 
lustrates the results. The DMC values with the JS wave 
functions are already quite good, because the cohesive 
energy, Ecoh, compares fairly well with the experimental 
values given in Ref. [2^, and because JS DMC calcula- 
tions are known to retrieve better than 90% of the cor- 
relation energy jsoj. Also, in both cases we observe the 
expected increase in cohesive energy with cluster size, 
describing the tendency for these clusters to grow [2^. 
However, the inset of Fig. [T] shows that the difference 
AE = E^^JfP - E-^J-"'", i.e. the nodal structure cor- 
rections, are about 30 meV/atom. 

The bonding properties of Li can be described quite 
accurately by replacing the core electrons with a pseu- 
dopotential, because the Is core states do not contribute 
significantly to bonding. For example, pseudopotential 
DMC calculations successfully predict the small binding 
energy for the LiPs molecule [3l| , in accordance with the 
results of more sophisticated DMC all-electron calcula- 
tions [13] ■ Another advantage of the pseudopotential is 
the possibility of relatively straightforwardly disentan- 
gling the valence characteristics of the A = Ay matrix 
in Eq. [T] from core contributions. This submatrix will 
also be seen below to help identify an RVB signature in 
the wave function. 

We consider the JS-LDA pseudopotential wave func- 
tion for Li clusters containing 2, 4, 8, and 20 atoms, in 
order to assess the impact of these one-body orbitals on 
the nodal structure. The inner Is core electrons were re- 
placed by the norm-conserving pseudopotential provided 



N 


pJS-LDA 
^VMC 

(Hartree) 


rpJAGP 
^VMC 

(Hartree) 


„2 

(eV) 


2 

(eV) 


2 


-0.20561(4) 


-0.21509(3) 


0.2921(9) 


0.068(5) 


4 


-0.20899(4) 


-0.22259(2) 


0.392(3) 


0.0477(9) 


8 


-0.21616(2) 


-0.231200(6) 


0.388(3) 


0.0662(5) 


20 


-0.21834(2) 


-0.2371(1) 


0.602(2) 


0.011(1) 



1^ 



TABLE III: VMC calculations of the JS-LDA and the JAGP 
wave function for Li clusters. The first two columns compare 
the total energy/atom. The variational JAGP results for Li2 
and Li4 are lower than the diffusion JS-LDA results from Ta- 
ble [III indicating superior performance of the wave function 
even at the variational level for the JAGP than can be ob- 
tained at the diffusion level for the JS-LDA. The final two 
columns show the variance of the local energy/atom. The 
JAGP shows a smaller variance than the JS-LDA. 



by Burkatzki et al. fss'l. The wave function was con- 
structed using Kohn-Sliam LDA orbitals obtained with 
the P Wscf code 3J] . The LDA calculations used a cubic 
simulation cell of 40 a.u. sides, a plane-wave cut-ofF of 70 
Ry, and included optimized Ji and J2 terms. A modified 
version of the DMC, the Lattice Regularalized Diffusion 
Monte Carlo (LRDMC) [35*1 method, was used, which 
allows the inclusion of a non-local pseudopotential in a 
consistent variational scheme. The JS-LDA DMC results 
are summarized in Table HIl In this case also, the JACP 
yields the lowest total energies. These results confirm 
that in the case of Li clusters the RVB nodal structure 
provides a correction to the cohesive energy, which tends 
to be more than 20 meV/atom. 

The signature of an RVB state can be directly iden- 
tified by analyzing the eigenvalues and eigenvectors of 
the matrix A = Xij in Eq. [1] Given the A matrix, 



we solve the eigenvalue problem ASu = fiu, where S is 
the overlap matrix between the local basis set orbitals, 
i.e. Sij — {4>i(j)j), and u are the eigenvectors representing 
the basis of natural orbitals, with corresponding eigen- 
values providing the geminal coefficients associated with 
the occupation numbers [36]. Anderson has shown [9] 
that a signature of the RVB state is given by a change 
of sign of the eigenvalues /i when passing from occupied 
to virtual states with a significant weight for the virtual 
states. This trend also helps to reduce double occupancy 
at lattice sites, as in the Heitler London limit. 

This diagonalization was performed for Li4, Lig and 
Li2o with pseudopotential. Fig. [2] shows the eigenvalues. 
Positive eigenvalues are in the top portion of the figure, 
while the negative eigenvalues are shown in the bottom 
on an expanded scale. The presence of the RVB signa- 
ture indicates a departure from the standard JS nodal 
structure, and by extension into the bulk, a departure 
from the Fermi liquid picture This departure can 



be explained as follows. If the generating geminal yields 
a standard JS ansatz, the occupied orbitals would have 
equal weights, and the remaining orbitals would not con- 
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FIG. 2; (Color online) Eigenvalues of the diagonalized A ma- 
trix, representing the occupation of the natural geminal or- 
bitals. Top half: eigenvalues for the first natural orbitals. 
Bottom half: expanded view of the eigenvalues for the virtual 
orbitals. Unlike a standard Jastrow-Slater wave function, the 
JAGP allows occupancy of the higher-level orbitals. A char- 
acteristic of the RVB is a sign flip when passing from the 
primary to the virtual orbitals. 



tribute. However, in the case of the RVB, the virtual 
orbitals have significant non-zero amplitudes, but with 
opposite sign. Contrary to the standard JS description, 
virtual orbitals, i.e. orbitals with index greater than A^/2, 
have a small but finite, negative eigenvalue, which is re- 
lated to the occupancy. 

The behavior of the eigenvalues in Fig. [2] confirms the 
fact that the AGP structure of the Li dimer studied by 
Blander et al. for various atomic separations [STj cap- 
tures general features which suggest a certain robustness 
of the RVB signature with respect to the exact geometry 
of Li clusters. 

It is interesting to note that the natural orbitals corre- 
sponding to the leading eigenvalues contain a significant 
p-character. An important 2p — 2p vr-contribution was 
observed in the AGP wave function for Li2 [37[, causing 
a departure from the standard JS nodal structure. Re- 
cent experiments find that Li impurities in an Al matrix 
produce an anomalous transfer from s- to p-like char- 
acter and thus constitute another example in which the 
standard Fermi liquid picture breaks down and proper- 
ties of the correlated inhomogeneous electron gas must 
be considered 



Table IIIII provides the results of Variational Monte 
Carlo (VMC) calculations, rather than DMC, of the JS- 
LDA and the JAGP wave functions, confirming the rel- 
atively high quality of the JAGP wave function at the 
variational level. A comparison of tables IIIII and [IT] in- 
dicates that our variational results using the JAGP are 
superior to the diffusion results using the JS-LDA for Li2 
and Li4. Table IIIII also shows the significantly smaller 
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variance of the local energy for the RVB wave function, 
indicating more efficient and less time-consuming calcu- 
lations [27], as well as a wave function that may be ap- 
proaching the exact solution. Interestingly, for the Li2 
molecule the JAGP is able to recover 95.7% and 99.45% 
of the correlation energy flG^ at the VMC and DMC 
levels, respectively, while the corresponding values with 
backflow corrections ^ are 87.79% and 97.1%. There- 
fore, in the present case, the RVB correlations appear to 
dominate other effects. Moreover, the possibility that the 
RVB Pauling structures provide an important contribu- 
tion in the configuration interaction (CI) expansion has 
also been suggested in a recent quantum chemical study 
for the neutral Li4 cluster [4o| . 

Finally, to directly measure nodal structure effects, we 
have utilized the limiting JS-VO wave function. In this 
way, we can estimate the energy contribution from oc- 
cupation of the virtual orbitals by forcing the pairing 
wave function to possess only the N/2 fully occupied or- 
bitals, with zero occupancy of the virtual orbitals. These 
calculations were performed for Li4, Lig and Li2o with 
pseudopotential, and in all three cases, VMC as well as 
DMC calculations find an improvement of at least 20 
meV/atom in the cohesive energy due to the RVB nodal 
structure, consistent with the DMC results of Table HIl 
As stated previously, the definition of JS-VO assures that 
the Jastrow factor, basis set, and functional form are 
identical in the JAGP and in the JS-VO wave functions, 
so that AEjivB = E;!^!^^ - E;!^-^° is a direct calcu- 
lation of the resonance energy. This confirms that the 
gain in energy is a consequence of the RVB orbitals, and 
not the result of a different Jastrow factor or a different 
single particle basis set. 

In conclusion, we have performed DMC calculations 
of Li clusters with an RVB guiding WF, utilizing a fun- 
damentally different nodal structure than the standard 
JS WF. We find that the RVB nodal structure is able 
to recover about 20 meV of cohesive energy /atom. Fur- 
thermore, we have identified a distinct RVB signature 
in an eigenvalue decomposition of the JAGP A-matrix, 
suggesting modifications of the electron occupation num- 
bers due to RVB effects. These results justify use of an 
approximation for AGP correlation effects in momentum 
density calculations 0] and might explain deviations from 
the Fermi liquid picture observed in recent high resolu- 
tion Compton scattering experiments on bulk Li. 
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